Graphene with geometrically induced vorticity 
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At half filling, the electronic structure of graphene can be modelled by a pair of free two- 
dimensional Dirac fermions. We explicitly demonstrate that in the presence of a geometrically 
induced gauge field, an everywhere-real Kekule modulation of the hopping matrix elements can 
correspond to a non-real Higgs field with non-trivial vorticity. This provides a natural setting for 
fractionally charged vortices with localized zero modes. For fuUerene-like molecules we employ the 
index theorem to demonstrate the existence of six low-lying states that do not depend strongly on 
the Kekule-induced mass gap. 
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Planar graphene and its geometrically related vari- 
ants offer a rich environment for exploring interesting 
physics [l| . The electronic properties of graphene are well 
modelled by a simple Hiickel model of nearest-neighbor 
hopping on a two-dimensional honeycomb lattice 0, Q ■ 
The low energy sector of the resulting band theory may 
be described by a pair of two-dimensional Dirac equa- 
tions. As a consequence, graphene is expected to ex- 
hibit phenomena more familiar in relativistic quantum 
theory, such as the Klein Paradox Recently, Hou 
et al. |5| demonstrated a relation between graphene and 
p-wave superconductors, where fractionally charged vor- 
tices can appear. Energetic considerations guided Jackiw 
and Pi @] to extend this model by inserting a gauge po- 
tential. Systematic study, including the demonstration 
of the fractional statistics of these vortices, was given 
in Nevertheless, the question of physically realizing 
these systems remained open. 

Spherical configurations of graphene are known as 
fuUerenes. The altered topology requires defects where 
twelve of the regular carbon hexagons are replaced by 
pentagons. The resulting frustration and curvature are 
accounted for in the Dirac equation by introducing a chi- 
ral gauge field and a spin connection 
0, [3] ■ The gauge-field fiux (but not the spin-connection 
curvature) enters into an index theorem where it is re- 
sponsible for the six anomalously low-lying states seen 
in the spectrum of Ceo and related molecules [l^ \v^ . 
Chemists have long surmised that in fuUerenes not all 
the nearest-neighbor hopping elements are equal. In nu- 
merical calculations with the Hiickel model it is found 
that the molecules can lower their electronic energy by 
undergoing a small Peierls distortion — usually called 
in this context a Kekule distortion. This change in the 
bond lengths introduces a scalar "Higgs" field into the 
Dirac equation. It has recently been observed, however, 
that vortices in a complex-valued Higgs fi eld can bind 
zero-energy mid-gap modes 



Here we demonstrate that the Kekule distortion Higgs 
field is not a simple scalar field, but is a section of the 
non-trivial gauge bundle. This means that a real-valued 
modulation of the hopping strengths in a fullerene will 
give rise to a complex-valued Higgs field that automat- 
ically contains vortices similar to those of Abrikosov or 
Nielsen and Olesen. Thus, fuUerenes provide a physical 
setting where vortices with fractionalized charge appear 
naturally. The number of zero modes bound by these 
vortices is equal to the number of zero modes required 
by the index theorem 2l|, [2^ . Indeed, they are the same 



modes. 

We begin with the description of a flat sheet of 
graphene. In the electronic tight binding (Hiickel) ap- 
proximation Q electrons hop on a two-dimensional hon- 
eycomb lattice with lattice constant a. This bipartite 
lattice can be decomposed into two triangular Bravais 
sub-lattices A a and A^. The hamiltonian can be written 



H=-J2 + Strk)albr+s, + H.c. 



(1) 



tEAa fe=l 



where the vector r gives the position of the vertices of the 
A A lattice and the vectors s^, k = 1,2,3 connect each 
site of the A^ lattice to the three adjacent sites of Ab- 
The fermionic operators Or and b^+s^ annihilate electrons 
in the sub-lattices Aa and Ab, respectively. At half- 
filling, graphene possesses two independent Fermi points 



at K± 



± 



) where the positive and negative en- 



ergy bands touch in conical singularities. For momenta 
near these Fermi points we can replace the full hamilto- 
nian with a Dirac approximation 



H = J *t(r)iJ*(r)d2 



(2) 



where H is the corresponding one-particle Dirac hamil- 
tonian and the spinor is given by \l/(r)"^ = {ua, Uh, Va, fb) 
with the a, b indices correspond to the A^ and A^ sub- 
lattices and M, V, correspond to the two Fermi points. 
For the particular case of the Kekule distortion given by 
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(3) 
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the Dirac hamiltonian takes the form 0] 



H 



tJ* ■ p (71$ 

cri<i>* cr ■ p 



(4) 



where p = —id, and we have taken an overall constant 
to unity. Here, cr = ((Ti,cr2) are the two-by-two Pauli 
matrices. Transformations of the form (7^ eg) 1 act on the 
Fermi point components (it, v) , while transformations of 
the form 1 ® aj act on the sublattice indices (a, b). Here 
we take $ to be initially real everywhere. 

We wish to compactify the Kekule distorted graphene 
sheet into a surface with the topology of a sphere. To 
do this it helps to make a change of basis. We exchange 
the role of A and B at the K+ Fermi point and rotate 
the reference frame at the K_ point by tt angle so that 
it coincides with the frame at K+. These two transfor- 
mations are effected by conjugating with the SU(4) ma- 
trix iJi ©as. To introduce curvature, we observe that the 
Kekule distortion leaves every third hexagon with no dou- 
ble bonds. We select one of these hexagons and, starting 
from its geometrical center, excise a wedge of opening 
angle 7r/3, the cuts passing through the centers of two 
lines of bonds. We then reconnect the dangling bonds 
to form a seam. This operation leaves us with a conical 
curvature singularity centered on a pentagon, but docs 
not cause a dislocation in the pattern of double bonds. It 
does, however, introduce frustration in the electron wave 
functions. This is because the identity of the A and B 
lattices is interchanged across the seam, and because an 
electron with wavefunction located near one Dirac point 
sees itself across the seam as a wavefunction belonging 
to the other Dirac point. To make the spinors continu- 
ous across the cut we therefore introduce a gauge-twist 
transformation 



* ^ [/*, with U = exp 



(ea -I- q) ■ dr (5) 



when the spinor is transported around the apex of the 
cone. Here a = A(ct2 ^ 1) is the non-abelian field with 
circulation around the pentagonal plaquette e f a - dr = 
7r/2((T2 <E) 1) and q is the spin connection with circula- 
tion ^q ■ dr = — 7r/6(l (g) (73). Taking into account that 
WdJJ = —ieAi{<T2®'^) — iqi we have that the momentum 
operator becomes p — i(V — leA) with V = 5 — iq. 
We can diagonalize the gauge field and simplify the Dirac 
matrices by the rotation {^q-^<^i^/^ (g) © CT3) giving fi- 
nally 



Ha 



-ie^^ai{V - ieA) <I> 

ie'^ai{V + ieA) 



(6) 



where the zweibeins, e], are introduced to make the 
Hamiltonian covariant in the induced curved space with 



metric 



where rjij is the flat metric. Note 



that, due to the geometric distortion, the initially real <!> 
has now acquired a phase factor of the form exp(i2e J A). 



Hamiltonian ([6]) corresponds to a spinor field defined 
on a curved manifold coupled to a chiral gauge potential. 
We observe that the gauge invariant Dirac operator. Ha, 
that we have derived coincides with the one recently writ- 
ten down by Jackiw and Pi ^ . In our case the pentagon 
is threaded by a quarter unit of gauge flux, and that the 
phase of the scalar field $ winds through only tt as we 
circle the defect. Each pentagon is therefore a rather 
singular "half- vortex" . To make the Higgs field smooth, 
we would need to join pairs of vortices and switch off 
the Kekule distortion along a line joining them (see Fig- 
ure [Ua,b)). We will discuss the significance of this issue 
later. 

We constructed our hamiltonian (jH) for the case of a 
single conical deformation. In order to make a surface 
with the topology of a sphere we formally need to in- 
troduce twelve pentagonal defects. Can we extend our 
surgical construction globally? It is known that we ob- 
tain a consistent assignment of double and single bonds 
composing a dislocation-free Kekule structure (a Fries 
structure) if and only if the fuUerene belongs to the fam- 
ily Cgo+efe, k = 0, 1, . . . that is obtained from a parent 
C2o+2fe molecule by "leapfrog" inflation [i^. This is ex- 
actly the family of molecules we obtain by the previous 
process, where each of the twelve pentagonal deforma- 
tions is surrounded by hexagons each possessing three 
double bonds. The total gauge flux through the twelve 
pentagons is 3 • 27r, and the total winding number of the 
Higgs field is six, composed of twelve half- vortices. This 
can be equivalently described by a monopole sitting in- 
side the molecule providing the net vortex flux [1, Q . 

At this point it is convenient to study the transfor- 
mation properties of the one-particle Dirac hamiltonian 
Ha- First, note that the different chiralities in the Dirac 
hamiltonian correspond to normalized solutions, which 
have their support on the different triangular sublattices 
of the honeycomb lattice. By comparison with H we 
deduce that the Dirac hamiltonian Ha acts on the vec- 
tor (lib, Mq, Ua, Wb). For convenience we change basis so 
that the spinor is given by (ua, I'a, Uh, Ufc). Then one 
has V^HA^b = —Ha where Fs = era © 1. This trans- 
formation maps positive energy eigenvectors of Ha to 
negative ones, i.e. r5\E'£;(r) = ^'_£;(r), while zero en- 
ergy modes are left unpaired. This property is known 
as sublattice symmetry and it gives rise to charge frac- 
tionalization and the emergence of zero modes when the 
scalar field has non-trivial vorticity [5|. Furthermore, by 
defining the operator Ei = 1 (g) cti we can show that 
TiiHa^i = H\. As a consequence [IS] the eigenf unc- 
tions of Ha satisfy Si^* — 5"* and in spinor components 
u* = Va and — Vh. This symmetry is known as the 
time reversal symmetry and it is present due to the choice 
of real hopping elements in hamiltonian ([1]) . 

We wish to apply the index theorem of Jackiw and 
Rossi [2l[, and E. Weinberg [S^l to our problem. This 
theorem was first introduced to study zero modes in a 
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relativistic analogue of superconducting vortices which is 
exactly the case here. As we shall see the number of zero 
modes provided from the index theorem can be computed 
by either counting the number of gauge field flux- units or 
by counting the net winding number of the Higgs field — 
the two numbers necessarily being equal because of the 
topology of the gauge bundle. 

The sublattice symmetry dictates that the non-zero en- 
ergy eigenfunctions come in ±E pairs, whilst zero-energy 
eigenf unctions (zero modes) can be chosen to be eigen- 
vectors of Tq. Thus, we can arrange the zero modes in 
terms of their chirality. Suppose that there are n+ zero 
modes with eigenvalue -1-1 and n_ with eigenvalue —1. 



dot 



Weinberg shows [22| that the Index(_ffDiiacj — 
is given by 



Index(i?Dirac 

e 



1 



F(Pr + (b dr ■ 

— f dr ■ dAig $. 
27r Jon 



$*D$ - $D$* 



(7) 



Here 



D$= (9-2ieA)$, D$* = (9 + 2ieA)$*, 

and the dfl integral is to be taken over a contour sur- 
rounding all the vortices. It is possible to use the consis- 
tency equation, {d — 2ieA)(f> = 0, of the scalar field that 
follows from Ai{r) = di(j){r). Then, it is easy to see that 
the index can be written in terms of the gauge field 



Index(_ffDii 



e 

71" Jn 



Fd^r. 



(8) 



From this form of the index of Huirac we deduce that 
one can take the scalar field continuously to zero, $ ^ 0, 
without changing the number of zero modes. 

As derived by Weinberg, the index theorem applies 
to the case of open boundary conditions. To apply the 
theorem to a spherical fuUerene with hamiltonian Ha re- 
quires the following step. We recognize that there is a 
net gauge-field flux through the sphere, and so the gauge 
and associated Higgs bundle are non-trivially twisted. 
We must therefore introduce two hemispherical patches 
and sew them together with a gauge-transformation 
that identifies d — leA. in the upper hemisphere with 
e^^^{d — ieA)e~^^^ in the lower, and similarly identi- 
fies $ in the upper hemisphere with e^"^''^$ in the lower. 
The phase ex in the U(l) group element e^'^^ will wind 
through — Gtt as we encircle the common bound- 

ary of the hemispheres. It is the twisting of the Higgs field 
that allows us to have a net number (six) of Higgs field 
vortices. If the field $ were an ordinary scalar instead of a 
section of a twisted bundle the net vortex number would 
necessarily be zero. Thus, the index theorem demon- 
strates that the mass term appearing due to the Kekule 
distortion does not actually destroy the zero modes. 



We have seen that the introduction of curvature 
through pentagonal defects automatically turns Kekule 
distorted planar graphene into a gauge theory of the form 
introduced by Jackiw and Pi. The defects form half- 
vortices in the Higgs field corresponding to the Kekule 
distortion. We can apply the Jackiw-Rossi- Weinberg in- 
dex theorem to find a lower bound on the number of low 
energy states on the curved surface. When applied to the 
leapfrog fuUerenes Cgo+efc we find that there should be 
six low-lying modes that are insensitive to the magnitude 
of the Kekule distortion. The singular nature of the vor- 
tices is a problem however, and numerical investigation of 
the leapfrogs C^n QO with a uniform distortion shows six 
low lying modes that do depend on the distortion (Fig- 
ures Ufa) and[3fa)). This is not surprising as in this form 
the Higgs field is discontinuous due to its half vorticity at 
each pentagon. To rectify it one must allow the vortices 
to pair up by joining them with cuts, that is regions with 
$ = 0. Moreover, the energy contribution from the scalar 
field <1> is given by / ii^2r{|D$p + y($*(f>)}, where F is a 
function with a minimum at <I>o with V'{^}^^a) — that 
enforces ^ to acquire mass. To avoid divergences of the 
energy the scalar field Q has to vanish polynomially as 
r — > 0. This forces us to consider vortices with enlarged 
size. As we shall see in the following these modifications 
in the configuration of the Kekule distortion of a fullerene 
molecule provide a spectrum that is in agreement with 
the prediction of the index theorem. 

Figures [2] and [3] depict the spectrum of C60 and C180 
molecules with Kekule distortion for the various cases 
of uniform and smoothed Higgs fields. The coupling of 
the single bond has been set to 1 and that of the double 
bonds is denoted by h. The case of uniform Kekule dis- 
tortion, presented in Figure [Ija), has been modified in 
order to ensure smoothness of the effective Higgs field. 
For that we first introduced cuts among the vortices, i.e. 
we removed double bonds along paths that connect pairs 
of vertices, as seem in Figure [Ijb). Moreover, we en- 
larged the size of the vortices by removing double bonds 
from the links attached to the pentagons, as seen in Fig- 
ure [Ijc). When the last process is performed on the C60 
molecule it removes all double bonds due to its small size. 

As expected from the index theorem there are six 
modes with near zero energy when the coupling h is var- 
ied. While for the case of uniform Kekule distortion these 
modes seem to be sensitive to variations of h (see Fig- 
ures Ufa) and [3Ja)) this is rectified by introducing cuts 
and by enlarging the vortices. These processes make the 
Higgs field continuous and the energy of the low-lying 
modes becomes even closer to zero with an energy that 
is insensitive to variations of h. Indeed, in Figures [2fc) 
and [3fc) there are six low-lying modes which are to a 
good degree insensitive to variations of h. The difference 
of their energy from being exactly zero is due to the small 
size of the system and it is expected to converge to zero 
when larger molecules are employed Q. 
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FIG. 1: The coupling configuration of the C60 molecule, 
where vortices reside on the pentagons, (a) The Kekule dis- 
tortion, (b) Cuts between vortices are introduced by replacing 
a double bond with single ones, (c) An enlargement of the 
vortices is introduced by removing all double bonds connected 
to the pentagons. For C60 this removes all double bonds. 




FIG. 2: The spectrum of the C60 molecule as a function of 
the double bond coupling, h, with single bond coupling equal 
to 1. There are six modes (two triplets) that are near zero 
energy, (a) With Kekule distortion, (b) With cuts between 
pairs of vortices, (c) With enlarged vortices, where all double 
bonds are removed. 




FIG. 3: The spectrum of the C180 molecule, (a) With Kekule 
distortion, (b) With cuts between pairs of vortices, (c) With 
enlarged vortices. 



To summarize, we explicitly demonstrated that the ef- 
fective gauge field induced in a graphene sheet, when it 
is geometrically deformed, couples non-trivially to the 
Higgs field induced from a distortion in the tunnelling 
couplings. As a result, the fullerene-like molecules have 
six half vortices that are fractionally charged [3| giving 
a natural setting where this topological effect appears. 
By employing the index theorem we demonstrated that 
there should be one zero mode for each pair of such vor- 
tices. While the effective gauge field is responsible for the 
presence of the zero modes the scalar field assures that 



the corresponding degenerate states are not locally dis- 
tinguishable implying their topological robustness (topo- 
logical degeneracy). In the case of fullerene-like molecules 
a study of the corresponding low- lying states showed that 
a Kekule distortion is energetically favorable actually re- 
ducing the overall energy of the electrons. 
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